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ABSTRACT 

Close binary systems undergoing mass transfer or common envelope inter- 
actions can account for the morphological properties of some planetary nebu- 
lae. The search for close binary companions in planetary nebulae is hindered 
by the difficulty of detecting cool, late-type, main sequence companions in bi- 
nary systems with hot pre- white dwarf primaries. However, models of binary PN 
progenitor systems predict that mass accretion or tidal interactions can induce 
rapid rotation in the companion, leading to X-ray-emitting coronae. To test such 
models, we have searched for, and detected, X-ray emission from three binary 
central stars within planetary nebulae: the post-common envelope close binaries 
in HFG 1 and DS 1 consisting of O-type subdwarfs with late-type, main sequence 
companions, and the binary system in LoTr 5 consisting of O-type subdwarf and 
rapidly rotating, late-type giant companion. The X-ray emission in each case 
is best characterized by spectral models consisting of two optically-thin thermal 
plasma components with characteristic temperatures of ~ 10 MK and 15-40 MK, 
and total X-ray luminosities ~ 10 30 erg s -1 . We consider the possible origin of 
the X-ray emission from these binary systems and conclude that the most likely 
origin is, in each case, a corona around the late-type companion, as predicted by 
models of interacting binaries. 
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Introduction 



Close binary systems have been fo und in 17 ± 5% of the central stars of planetary 
nebulae (CSPNe) (IMiszalski et al.l 120091 ). These systems are typically comprised of a hot, 
compact white dwarf (WD) or pre- WD primary and a low mass, late-type, mai n sequence 
companion and, hence, a re likely the products of common envelope interactions (IPaczynski 
19761 : llben fc Livid Il993l ). In such systems, the primary star overfills its Roche lobe as it 
evolves onto the red giant branch (RGB) or asymptotic giant branch (AGB), engulfing the 
main sequence companion and forming a common envelope (CE). Drag forces in the dense 
CE drain orbital energy from the binary, bringing the two stars closer together. The CE 
phase ends when either the companion spirals onto the primary's core; or the envelope is 
ejected. In the latter case, the residual is a close binary called a post-common envelope 
binary (PCEB). Although PCEBs in PNe are thought to derive from CE interactions that 
took place during the AGB phase of the primary or s econdary, other post-CE bi naries may 
have originated in CE interactions on the RGB (e.g., ISchreiber fc Gansickdl2003l ). 



The ejected CE can account for the morphological properties of some PNe. iBond &: Livio 



(119901 ) argue that PNe with dense equatorial waists, often called butterfly PNe, form from 
a CE ejection immediately following the primary's early AGB phase, whil e elliptical PNe 



form from a CE ejection foll owing a later AGB pha se. Recent observational ( IMiszalski et al. 



20091 ) and hydrodynamical (IRicker fc Taaml 120081 ) studies appear to bear out these asser- 
tions. Indeed, within the PNe community, there is momentum building behind th e idea that 
the e ntire diverse population of PNe may be the result of binary interactions (IDe Marco 
20091 ). 



In most studies of PCEBs emphasis is placed on the evolution of the primary star into 
a mass accreting WD, i.e., a cataclysmic variable (CV). However, in this paper we are more 
interested in the evolution of the late-type, secondary companions. The rapid rotation in 
youn g, single stars generates magnetic fields strong enough to power bright X-ray coro- 
nae ( jMathioudakis et al.lll9951); however, as thes e stars age, their magnetic activity decays 



see reviews m 



Glide! 



2004 



Giidel fc Nazel 120091 ). Hence, the old, late-type, companions in 



PCEB are not expected to exhibit strong coronal activity unless t heir rotation rates have 



increased by accretion of angular momentum during the CE phase (I Jeffries fc Stevens! 11996 



Soker fc Kastnerl 120021 ). Although the degree to which the companion m ay undergo such ac- 
cretion is uncertain (e.g.. lHjellming fc Taamlll99ll : IRicker &: Taamll2008l ). there is compelling 
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evidence from extre me UV and X-ray observa tions of PCEBs that such a process occurs. 
From EUV spectra, ISanz-Forcada et al.l (120031 ) find evidence for strong coronae around the 
cool companions in close binaries, while abundance anomalies determined from high resolu- 
tion X-ray spectra of the close binary V471 Tauri show that the K dwarf companion must 
have accrete d the highly processe d and enriched material liberated by the primary on its 
RGB ascent (brake fc Sarnall2003h . 



By analogy with main sequence stars, coronae around the cool companions in PCEBs 
should be characterized by thermal plasmas from a few MK to tens of MK and therefore 
should emit hard X-rays with energies > 0.5 keV. ISoker &: Kastnerl (120021 ) argue that for 



CSPNe, a plasma temperature above 10 MK and Lx > 5 x 10 29 ergs s _1 would serve as 
compelling evidence for a "reborn" corona around a spun- up companion. There is little or 
no coronal contribution expected from the hot primary, since coronal activity from hot WDs 
or pre- WDs with effective temperatur es >30 or >60 kK, respectively, is difficult t o maintain 



when photospheric convection ceases (jBohm fc Cassinellil Il97ll ; iGroth et al. 



1985). Further- 



more, if the CE is ejected, the magnetic field is ejected along with the CE (INordhaus et al. 
20071 ). Any remainin g magnetic field will s low down the rotation of the remnant AGB core 
via magnetic braking (IBlackman et al.ll200ll ). resulting in decreased coronal activity. Indeed, 
O'Dwyer et al.l ( 120031 ) discovered numerous sources of hard X-ray emission from WDs and 
argued that in all but two cases, this hard X-ray emission is a result of either accretion from 
a companion onto the WD or coronal activity on the companion itself. 

Thus, the potential exists to use X-ray observations as a tool for detecting otherwise un- 
seen binary companions to CSPNe. As an initial test, we targeted two binary systems, HFG 
1 and DS 1, with the Chandra X-ray Observatory (CXO), and analyzed archival data from 
the serendipitous observations by CXO and the XMM-Newton X-ray Observatory (XMM) 
of the binary CSPNe in LoTr 5. All three CSPNe were detected in these X-ray observa- 
tions. In this paper, we present analysis and interpretation of the X-ray emission from these 
three binary CSPNe. In Section 2, we summarize the properties of these binary systems; 
in Section 3 we describe the data and analyses; and in Section 4, we discuss the possible 
interpretations of the X-ray emission and their implications for models of binary interaction 
and PN shaping. 



2. Target CSPNe 

The binary systems we consider are the PCEBs in the PNe HFG 1 and DS 1, and the 
evolved binary in LoTr 5. Throughout the text, we will generally refer to both the central 
star and PN by the PN name. In each case, the primary, and presumed origin of the CE 
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and PN, is a hot (T e ff > 60 kK) O-type subdwarf (sdO) evolving towards a WD, while 
the secondaries are cool, late- type stars (see Tabled]). These systems are the precursors to 
cataclysmic variables (CV), but do not exhibit characteristics of active accretion, and they 
feature strong reflection effects, i.e., irradiation of the secondary by the primary. 



2.1. HFG 1 and DS 1 



Properties of the PCEB syste ms in the PNe HFG 1 (V664 Cas) and DS 1 (KV Vel, 



LSS 2018) have been summarized by lDe Marco et al.l (120081 ) and are listed in Table [TJ Here, 
we briefly highlight the key pro perties. The binary orbital periods are well-determined due 



to strong reflection effects (e.g. lExter et al. 



2005 



Hildi tch et al.l Il996l ) . The binary period 



of V664 Cas (within HFG 1 ) is 14 hours (IGrauer et al.l 119871 ) . Assuming a primary mass 
~ O.6M , lExter et al.l ( 120051 ) determined a r ange of acceptable fits f or the secondary mass 
(M 2 ~ 0.4 — I.IMq). The spectral analysis by IShimanskii et al.l (120041 ) suggests the chemical 
composition of the secondary in V664 Cas is near solar, but with some anomalies (overabun- 
dances of nitrogen, magnesium, and silicon) that are interpreted as enriched prim ary material 



accret ed by the secondary during the CE phase. Distance estimates to HFG 1 by lExter et al. 



d2005h range from . 31 to 0.95 kpc, which is consistent with the statistical determination 
by iHeckathorn et alJ ( 19821) of .4 kpc and a distance of 0.6 kpc found from the PN surface 
brightness relation ( IFrew II2008I ). Here, we adopt a value of 0.6 kpc. 



ditch et al. 



The binary period of KV Vel (within DS 1) is 8.5 hours (jDrillinglll985l ). iHi 
1990) determine M 1 = O.63M and M 2 = O.23M for the components of KV Vel. 
19841 ) used the color excess measured from KV Vel to estimate a distance range of 0.5 to 



Schonberner fc Drilling 



1.4 kpc. We ad opt the distance of 0.7 kpc obtained via the PN surface brightness relation 
bv lFrewl (b()08h . 



The PNe HFG 1 and DS 1 are large and faint (Bond fc Liviclll990). HFG 1 shows strong 



evide nce o f interacting w i th the interstellar medium ( IHeckathorn et al 



1982: Heckathorn Fesen 



19851 ) and lBoumis et al.l ( 120091 ) recently discovered a trail of shocked material behind the star 
as it plows throu gh its local environm ent. The PN DS 1 suggest s a late elliptical PN with 
blown out edges (IBond fc Livid 1 1 9901 ). and iMiszalski et al.l (120091 ) point out the presence of 
low ionization structures that suggest an outflow along the E-W axis of the PN. 
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2.2. LoTr 5 



The binary system in the PN LoTr 5 (IN Co m, HP 112313) is comprised of an sdO pri- 
mary a nd a rapidly rotat i ng; G5 III-IV companion jFeibelman fe Kaler![l983l ; IStrassmeier et al. 



19971 ). IJasniewicz et a 



stars (IBopp fe Stencel 



( 19961) considered the similarities with rapidly rotating FK Comae 



19811 ) . arguing that FK Comae systems end the CE phase as a coa- 



lesced star, whereas LoTr 5 ended the C E phase as a close binar y after ejecting its envelope. 
As the binary orbital period is unknown, IJasniewicz et al.l ( 119961 ) consider a few alternatives: 
a) a very short period of ~ 1 day, b) a moderately short period of a few days, and c) a wide 
binary with a period of a few years. The radius of the giant star (~ V2R Q ) rules out the close 
separation required for a very short period binary, but beyond this lower bound, there is 
little constraint on the binary orbital period. The observed photometric variation of 5.9 days 
is attributed to the rotation of t he giant companion (see detailed discussion on the history 
of orbital solutions for LoTr 5 in IStrassmeier et al.l 119971 ) . This interpretation is supported 
by the presence of Ca II H & K emission lines that are indicative of chromospheric activ- 
i ty due to rotation, and which imply a projected rotational velocity of v sim = 67 km s -1 
(IStrassmeier et al.lll997| ). 



Long sl it spectra and narrow band images of the large (~ 6(5 x 1(7), faint PN are 
modeled by iGraham et al.l (120041 ) as a bipolar nebula inclined 17° to the line of sight. If 
the nebula is formed from the ejection of a CE , then the orbital plane is likely to coincide 



with the bipolar axis (e.g. iBond fe Livid Il990l ). Such a low binary inclination would then 



explain the inabil i ty to observe radial velocity variability from the binary system. However, 
Strassmeier et al.l (119971 ) found that the inclination of the rapidly rotating giant is i ~ 45°. If 
one takes the orbital inclination to be identical to that of the bipolar lobes, t hen this suggests 



that t he binary and its components are not coplanar. This is considered in IJasniewicz et al. 



(119961 ) as evidence for a longer orbital period , since the components are not close enough to 
become synchronized via tidal interactions. iThevenin fe Jasniewica (119971 ) argue that the 
high projected equatorial velocity suggests LoTr 5 is a wide binary system, in which the 
rapid rotation inferred for the giant companion is due to accretion from the AGB wind and 
does not require a CE phase. 



Distance estimates to LoTr 5 are widel y discrepant, ranging fro m 0.5 to 6.9 kpc (IGraham et al 



20041 ); the clo ser dis t ance is favored (e.g. IStrassmeier et al.l 1 199 71 ). We adopt a distance of 
0.5 kpc, after iFrew I (120081 ). The high galactic latitude (+88?46) of this PN renders inter- 
stellar reddening negligible. 
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3. Data and Analysis 
3.1. Observations 

We acquired targeted observations of HFG 1 and DS 1 on the back-illuminated S3 chip 
of the ACIS detector array onboard CXO. A summary of the CXO observations of HFG 1 
and DS 1 is presented in Table |2j In these two on-axis CXO observations we detect X-ray 
point sources centered on the positions of the central stars in HFG 1 and DS 1. There is 
no evidence for diffuse X-ray emission from either nebula (see also Section 4). These CXO 
observations were prepared and analyzed according to d ata analysis threads a ccompanying; 



the Chandra X-ray Center CIAO software (version 4.1) (IFruscione et al.l 120061 ). We find no 



high-background periods during the CXO observations of HFG 1 and DS 1. 

Archival data of the serendipitous observations of LoTr 5 from both CXO and XMM were 
obtained through the High Energy Astrophysics Science Archive Research Center, a service 
of the Astrophysics Science Division at NASA's Goddard Space Flight Center and the High 
Energy Astrophysics Division of the Smithsonian Astrophysical Observatory. A summary 
of the XMM and CXO observations of LoTr 5 is presented in Table [2J LoTr 5 lies 7f3 off 
axis in the XMM observation (ObsID 0012850201; 2002 June 06) of the galaxy group NSCS 
J125606+255746 and is detected on all three European Photon Imaging Camera (EPIC) 
detector arrays (pn, MOS1, and MOS2). The XMM observations were operated in Full- 
Frame Mode with the thin filter. Data were reprocessed completely from the Observation 
Data Files using the XMM- Newton Science Analysis Software (SAS) package version 7.1.0 
with the calibration files available in Current Calibration File Release 241 (XMM-CCF- 
REL-241). In the XMM observations of LoTr 5, the detected emission is consistent with 
a point source at the off-axis position of the CSPN. There is no evidence of diffuse X-ray 
emission associated with the nebula (see also Section 4). The Chandra observation (ObsID 
3212; 2002 December 04) of the same galaxy group field was obtained six months after 
the XMM-Newton observation. In this CXO observation, LoTr 5 lies ~ 8f6 off axis in the 
CXO observation, just within the field of view of the Advanced CCD Imaging Spectrometer 
(ACIS) detector array, on a CCD (the front-illuminated S4) that lies adjacent to the prime 
imaging CCD used to target the galaxy cluster (the back- illuminated S3). Since the ACIS-S 
chip array is optimally aligned for use with the curved focal plane Rowland circle for grating 
spectroscopy, the S4 image of LoTr 5 is very out of focus. Hence, it is difficult to ascertain 
the nature of the emission; however, its spectrum is useful. The CXO data obtained for LoTr 
5 were prepared and analyzed according to data analysis threads accompanying the CIAO 
software. We filtered the XMM observations of LoTr 5 for high background periods and bad 
events using the standard filters for the imaging mode observations. No high background 
periods were found during the CXO observation of LoTr 5. A summary of the observations 
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of LoTr 5 is presented in Table |2j 



3.2. Spectral Analysis 

The source and background spectra from the on-axis X-ray point sources in HFG 1 and 
DS 1 were extracted using the psextract task in CIAO. The source and background spectra 
from the off-axis, serendipitous XMM and CXO sources at the position of the central star of 
LoTr 5 were extracted from source regions determined by the 90% encircled energy radius 
for a point source at the appropriate source off-axis angles. For the XMM observations of 
LoTr 5 the 90% encircled energy radius is 45" and for the CXO observation the 90% encircled 
energy radius is 10". For each source spectrum, we generated target- and observation-specific 
response matrix (RMF) and ancillary response (ARF) files. This extraction of the source 
and background spectra and responses accounts for the broadening and vignetting due to 
the off- axis angle of the serendipitous observations of LoTr 5. 



All spectral fits were performed with XSPEC version 12.3. lx (lArnaudlll996f ). We use 
reduced xl statistics to distinguish best fit models, where v represents the number of degrees 
of freedom. When multiple models give acceptable xl values (xl < 2) we use goodness-of-fit 
tests (goodness in XSPEC), to determine the more appropriate model. This goodness-of-fit 
test performs Monte-Carlo simulations of the spectrum as drawn from the best-fit model, 
and calculates the xl statistic for each simulation. A value of ~50% is returned if the data 
is accurately described by the set of simulated spectra, while values closer to 100% suggest 
the data are not well described by the model. We performed 500 simulations per test. Due 
to the small number o f counts detected for HFG 1 and DS 1, we used Churazov weighting 



(IChurazov et al.l 119961 ) on the unbinned spectra. Such weighting is preferred for sources 
with low count rates and many empty spectral channels. Weighting is assigned with respect 
to the average of the surrounding channels. The best-fit spectral parameters (and their 
90% confidence levels) and the resulting inferred X-ray emission properties are presented in 
Table El We now discuss specific aspects of the spectral fitting for each source. 



3.2.1. X-ray Spectrum of HFG 1 

We attempted to fit single-temperature-component (1-T) and two-temperature-component 
(2-T) optically thin thermal plasma (mekal) models to the X-ray spectrum of HFG 1. At- 
tempts to determine the intervening absorption from the model fits failed to constrain the 
absorbing column density (Nh), so we adopted an intervening absorption, Nh = 2.4 x 
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10 cm , corresp onding; to the value of Eb- 



v 



0.5 mag estimated to the central star 



(jExter et al.ll2005l ). The absorbed 1-T model can only fit the data when the plasma abun- 
dance is allowed to vary. The best-fit 1-T model gives Y 2 n ~ 0.80 and requ ires an abundance 
of 0.06 times solar (with solar values given by lAnders fc Grevessdll989f ). with 90% confi- 
dence range of 0.02 t o 0.14 solar. However, this abundance result is inconsistent with the 
conclusions drawn by IShimanskii et al.l ( 120041 ) who found stellar abundances near or slightly 
above solar. The goodness-of-fit test yields 98% of the simulations with a better xl value, 
suggesting the flexibility introduced by freeing the model abundance parameter simply allows 
the model to fit features that are due to noise. Adding a second thermal plasma component 
improves the overall appearance of the fit and produces a similar xl value of xtio ~ 0.79. 
The goodness-of-fit test yields 68% of the simulations with a better xl value, however, indi- 
cating that the 2-T plasma model is a better description of the data than the 1-T, variable 
abundance model. The parameters of this best-fit 2-T model are presented in Table |3] and 
the model is overlaid on the source X-ray spectrum in Figure [TJ 



3.2.2. X-ray Spectrum of DS 1 

As in the case of HFG 1, the intervening absorption to DS 1 could not be constrained. 
We adopted Nh = 8.3 x 10 20 cm" 2 based on the va lue of Eb-v ~ 0.15 mag estimated 



toward the central star (ISchonberner fc Drilling! Il984j ) . The best-fit 1-T model with solar 
abundances gives X512 ~ 0-49 for a hot component at Tx ~ 14.3 MK, while the best-fit 
2-T model gives X510 ~ 0-49 for a similar hot component at Tx ~ 14.5 MK and a cooler 
component at Tx ~ 3 MK. The best-fit 1-T and 2-T models are indistinguishable, based 
solely on their xl values. The X-ray fluxes determined from the best-fit 1-T and 2-T models 
agree within their 90% confidence ranges. We find from visual inspection, however, that the 
2-T model locally improves the fit at the soft end of the spectrum. Hence, we adopt the 
2-T model as the best fit. Additionally, the goodness-of-fit tests for the two models shows 
a modest improvement, from 85% to 63%, when the second component is added. Adopting 
the 2-T model over the 1-T model does not effect the overall interpretation of the X-ray 
emission observed from the central star of DS 1. The parameters of the best-fit 2-T model 
are presented in Table |3j and the model is overlaid on the source X-ray spectrum in Figure |2j 



3.2.3. X-ray Spectra of LoTr 5 



We attempted to fit single-temperature models with wide ranges of absorption and 
elemental abundances to the XMM and CXO spectra of LoTr 5, but all such fits resulted 
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in x^3o ~ 3 — 6. The best-fit models of the X-ray spectra are instead comprised of two, 
optically thin, thermal plasma components. Attempts to add foreground absorption to the 
model did not appreciably improve the fit and did not constrain the absorbing column 
density. This is consistent with the high galactic latitude and proximity of LoTr 5, which 
suggest there is negligible intervening absorption. The parameters of the best fits of the 
X-ray spectra using 2-T models are presented in Table [3] and overlaid on the source X-ray 
spectra in Figures |3] & HI The difference between the apparent flux levels measured from the 
XMM and CXO observations is larger than t he cross-calibrat ion discrepancy between EPIC 
observations and CXO ACIS-S3 observations (jSnowdenl 120021 ). however, there are no specific 
cross-calibration data comparing XMM EPIC observations and CXO ACIS-S4, where the 
X-ray emission from LoTr 5 is detected. To facilitate direct comparison, we present the 
unfolded spectra from the CXO ACIS-S and XMM EPIC pn detectors together in Figure 
In the absence of any significant cross-calibration discrepancies, it appears that in the six 
month interval between the XMM and CXO observations the brightness increased across 
the entire energy range, indicating that the two spectral components arise from a common 
physical process. 



3.3. Timing Analysis 



To test for variability, we performed Kuiper tests (IPress et al.lll992f ) to assess whether 
the cumulative distribution of the barycenter-corrected source and background photon arrival 
times can be described by a Poisson process. These tests are independent of temporal bin size, 
but do not provide information on the period or amplitude of variability. Where the Kuiper 
test suggested deviation from a Poisson process, we constructed background-subtracted light 
curves using for multiple temporal bin widths. These background-subtracted light curves 
were extracted using the same source and background regions used for spectral extraction 
(using the ltd option of the CIAO tool dmextract). The low source count rates and relatively 
short temporal windows do not permit us to effectively use Fourier period detection methods 
or create phase-binned spectra. Hence, we restricted our assessments of source variability to 
the Kuiper test, comparision with the background, and visual inspection of the light curves. 



3.3.1. HFG 1 and DS 1 

The Kuiper test gives a 98% and 87% probability, respectively, that the photon arrival 
times of HFG 1 and DS 1 can be ascribed to a constant Poisson process. The minimum 
temporal bin widths to obtain a signal-to-noise ratio of 3 are 670 s and 3300 s for HFG 1 
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and DS 1, respectively. Below this signal-to-noise ratio the source count rate approaches the 
background rate and it becomes difficult to discern source variability from noise fluctuations. 
At these temporal bin widths, the rms noise is on par with that expected from a Poisson 
distribution. Improving the signal-to-noise ratios by widening the light curve bins to half the 
exposure times for HFG 1 and DS 1 confirms the Kuiper test result that HFG 1 produces a 
steady count rate, and suggests that DS 1 may exhibit variability (Figured]). 

3.3.2. LoTr 5 

For LoTr 5, Kuiper tests yield only 16% and 2% probabilities, respectively, that the 
X-ray sources detected in the XMM and CXO observations can be described by a constant 
Poisson process. In both observations, this deviation between the source and a simulated 
constant Poisson process is not present in the background. The minimum temporal bin 
width for a signal-to-noise ratio of 3 is 200 s and 320 s for the XMM and CXO observations, 
respectively. In the XMM observation, for temporal bin widths greater than 1000 s, the rms 
noise of the source is consistently 1.3 times that expected from a Poisson distribution. In 
the CXO observation, the rms noise of the source increases to 1.8 times that expected from 
a Poisson distribution. Since there are too few counts to construct phase-binned spectra, we 
cannot make proper inter-observatory comparisons of fluxed light curves. The background- 
subtracted light curves are presented in Figure H These light curves and the results of the 
Kuiper tests suggest that the X-ray source in LoTr 5 is modestly variable. 

4. Discussion 

As noted in Section 1, we anticipate that the X-ray emission from the central stars of 
DS 1, HFG 1, and LoTr 5 is due to the presence of coronal activity associated with the late- 
type companions in these likely or candidate PCEBs. However, there are other mechanisms 
which may give rise to thermal plasma emission in the range of temperatures we observe. 
We consider these alternative interpretations before discussing the coronal interpretation. 

4.1. X-ray Emission from the WD 



White dwarfs can be sources of X-ray emission, if they have low opacities or high effective 
temperatures (e.g. IJordan et al.lll994t iMarsh et al.Ml997t iMotch et al.lll993l ). However, WD 
photospheric emission is much softer than the X-ray emission detected from the three CSPNe 
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considered in this paper. Theoretically, cool WDs with connective envelopes may generate 



magnetic fields and coronae and emit X-rays (lSerberlll990l; iThomas et al.lll995l) . however, no 



such emission has been detected from field WDs ( jWeisskopf et al.l 120071 ) and the primaries 
in the three CSPNe considered here are too hot to sustain coronal activity. Hard X-ray 
emission from single WDs at Tx similar to those of the three CSPNe considered here is very 
rare; only two cases are known: KPD 0005+5106 and the central star of the Helix Nebula 
(jO'Dwyer et al.ll2003l ; iGuerrero et al.l 1200 ll ; IChu et al.ll2004l ). It is possible, but statistically 
unlikely, that the X-ray emission detected from one or all of the three CSPNe considered 
here is of the same (unknown) origin. 



4.2. X-rays from Colliding Wind Shocks 



Shocks from colliding winds associated with these CSPNe may come in two forms: 
(1) a fast stellar wind colliding with a slower, previously ejected, stellar wind — as in the 
inter acting stellar wind (ISW) theory that successfully explains the emissivity patterns of 
PNe (IKwok et al.lll978l ); and (2) collision of the two winds generated by the components of 
low-mass binary systems. Both types of energetic wind interactions can create X-ray emitting 
shocks with plasma temperatures in the range 1-8 MK and w ith X-ray luminosities from 10 31 



to 10 33 erg s 1 



e.g. 



Muerset et al.lll997l ; iKastner et al.ll2008l ). A significant contrast between 



the two types of colliding winds is the scale at which each occurs. The ISW collisions form 



large hot bubbles that fill the PN cavity on the order of a few thousand AU (e.g. IKastner et al. 



2000l ; iMontez et al.ll2005l ). while colliding-binary winds occur in a thin, conical interaction 
region with a vertex loca ted at the stagnation point of the two winds, typically on the order 
of the binary separation (iMuerset et al.lll995l ). 



In the targeted, on- axis, CXO observations of HFG 1 and DS 1, and the off-axis XMM 
observation of LoTr 5, there is no evidence for any extended X-ray emission. In the case of 
NGC 40, the X-ray emission is of such a low surface brightness that it was only rev ealed 
via blind spectral extraction from the expected emission region (IMontez et al.l 120051 ) . We 
performed similar blind spectral extractions for the X-ray observations of the large PNe in 
LoTr 5, HFG 1, and DS 1 and found no evidence for diffuse, low surface brightness X-ray 
emission. Furthermore, the X-ray luminosities determined from the X-ray spectral fits to the 
point source emission in these three, large PNe (Tabled]) are fainter than the faintest diffuse 
X-ray emitting plas ma in PNe (~ 10 31 erg s _1 ) and the plasma temperatures are hotter (see 



Kastner et al 



20081 ) 



The emitting volumes of the unresolved X-ray emission can be explored through the 
emission measures determined from the X-ray spectral fits. The emission measure (EM = 



J n e nndV ~ n 2 e V) of large, diffuse, X- ray emitting plasmas in PNe suggest plasma densities 
of a few 100 cm -3 (jSteffen et al.ll2008f ). where as the emission meas ure from colliding binary 
winds suggest densities of 1 — 5 x 10 6 cm~ 3 (IMuerset et al.lll995l ). For a typical emitting 
volume characteristic of colliding binary winds (~ 10 41 cm 3 ) the emission measure determined 
from the X-ray spectral fits (Table |3]) suggest the density is ~ 10 6 c m" 3 . This density i s 



similar to that inferred for colliding binary winds in symbiotic stars (IMuerset et al.l 119971 ). 
However, the X-ray luminosities reported in Table [3] are one to two orders of magnitudes 
fainter than such systems. Additionally, in HFG 1 and DS 1, there is no ev idence for a wind 
fast en ough (> 500 km s _1 ) to account for colliding-binary wind shocks. iModigliani et al. 
( 119931 ) reported a fast wind speed of 3300 km s _1 from their analysis of the IUE high 
resolution UV spectrum of LoTr 5 that potentially could explain the temperatures determined 
from the X-ray spectra of LoTr 5, although the evidence for such a fast wind is disputed 
(Guerrero; private communication). 



4.3. X-rays from Accretion-related Processes 



Although there is no evidence for accretion disks in prior studies of HFG 1 and DS 1, 
given that these systems are candidate pre-cataclysmic variables, we should consider whether 
the X-rays from the core of each PN could be due to a h ot accretion disk or from accretion 
onto the compact sdO primary. From lFrank et al.l ( 120021 ). the maximum temperature for an 
accretion disk around a compact object (in this white dwarf) is given by 



T 



120kK 




M 



R 



1 x 1O- 7 M year 



-i 



O.Oli?, 



1/4 



Considering the primary star masses and radii listed in Table (U the maximum disk tem- 
perature for these three binary systems would exist around the compact primary in LoTr 
5. However, this temperature, T max ~ 40 kK, is much too cool to account for the detected 
X-ray emission, hence, the X-ray emission cannot arise from a hot accretion disk, but this 
does not preclude the existence of a disk, nor a disk as a source of accretion material. 

If material is falling onto the compact object from Roche lobe overflow (RLO) from the 
companion, or from the inner region of an accretion disk. Such material can reach velocities 
up to the free fall velocity, Vff = (GM /R) 1 ^ 2 , forming shocks in the boundary layer (BL) 
between the flow and compact object. From iFrank et al.l (120021 ). the temperature of the 
shocked material is given by 



13 



The resulting implied temperature for the compact companion in LoTr 5 would heTff ~ 90 
MK. This is higher than the temperature determined from the X-ray spectral fits, though 
not high enough to confidently dismiss as the origin of the X-ray emission. Meanwhile, 
the implied temperature for the compact companions in HFG 1 and DS 1 is in the range 
Tff ~ 8 — 10 MK, which is consistent with the range we determine from the X-ray spectral 
fits. 

Since an accretion disk and/or inf ailing material can imply the presence of RLO, we 
consider the possibility of RLO in the three binary CSPNe. We estimate the Roche lobe 
(Rl) radius using the approximation given in lEggletonl (119831 ). 



Rl 
A 



0.49g 2 / 3 



0.6g 2 / 3 + ln(l + gV3) 



(3) 



where A is the orb ital semi-ma j or ax is and q = M1/M2. We find that a few of the possible 
solutions given by lExter et al.l (120051 ) for the binary system parameters in HFG 1 suggest 
that the main sequence companion can fill its Roche lobe, but the solution provided by 



Shimanskii et al 



2004 J) sugg e sts it cannot (see Table [T]). The binary system parameters 



determined by iHilditch et al.l (119961 ) suggest that the main sequence companion in DS 1 is 
unable to fill its Roche lobe and, therefore, is not capable of actively supplying material for 
accretion streams or accretion disks. The lack of support for RLO suggests that accretion is 
likely not the source of the observed X-ray emission from HFG 1 and DS 1 . 

For LoTr 5, if we assume synchronized orbital and rotational periods (P or b = P ro t), 
then the binary separation, A, is ~ 16-R© and the Rl ~ In this scenario, ey en at 

the lower limit on the giant's radius (~7.8 R Q according to IStrassmeier et al.lll997l ). the 
giant can lose mass to the sdO primary via RLO. If the orbital period of LoTr 5 is longer 
than 13 days, the giant is unlikely to be undergoing RLO. H owever, there is eviden ce of 
short-term variability (perhaps rapid flickering) in the optical (IStrassmeier et al.lll997l ) and 
longer-timescale variability in the X-ray (Figure [7]). B oth forms of variability could be 
indicative of accretion on the sdO star, a la CV system s (lBruchlll992h. A period of ~ 0.25 



days was reported by lKuczawska fc Mikolajewskil (119931 ) for LoTr 5. IStrassmeier et al.l (119971 ) 
searched for a similar periodicity but could not find one, leading them to suggest the period 
found by iKuczawska fc Mikolajewskil (119931 ) was an artifact or a time variable phenomenon. 
The origin of t his putative 0.25 day pe riod could be a hot spot on an accretion disk around 
the sdO, or, as IStrassmeier et al.l (119971 ) suggest, a reflection effect off an as- yet undiscovered 
third component in Lo Tr 5. The X-ray lightcurves i n Fig ure [7J do not clearly support the 
periodicity proposed by IKuczawska &: Mikolajewskil (119931 ). but we cannot be certain from 
our data. Further observations are required to determine if such a periodicity can be found 
in the X-ray source in LoTr 5. 
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4.4. X-rays from Internal O-star Wind Shocks 



The shocked winds of and B stars display log Lx/L-^^ ~ —7 (jCassinelli et al.lll994l . 

198fl ). If we attribute the detected X-ray emission to the primary in HFG 1, t h en lo g 



Lx/L-^ Y ~ —6.1 ± 0.2, which is higher than that found in ICassinelli et all ( 119941 . Il98ll ). 
For the primary in HFG 1, a shock velocity of ~ 1200 km s -1 relative to the wind veloc- 
i ty is required to p roduce the measured 20 MK plasma. From a radiative driven winds 
(jOwocki et al.lll988l ). this would suggest the wind speed is at least twice the shock velocity, 
i.e., > 2000 km s _1 . As of yet, no such wind has been found in HFG 1; hence, it is unlikely 
that the X-ray emission arises from shocks in the wind from the primary. Alternatively, 
the observed value of L x / -^ do i could suggest a magnetically active primary, but as we have 
stated, such a scenario is difficult to reconcile in a star with such a high effective temperature. 
Similarly, attributing the X-ray emission from the point source in DS 1 to the primary leads 
to log Lx/^boi ~ —6.6 ± 0.3, again suggesting a fast shock velocity (~ 1000 km s -1 ) or a 
high level of magnetic activity in a star that is unlikely to support active magnetic fields. 



4.5. X-rays from Coronal Activity Associated with the CSPNe Companions 



Coronal X-ray lumino sity is strongly corre lated with rotation for late- type main se- 
quence stars (see review by lGiidel fc Nazell2009l ). This correlation is believed to arise, as in 
our sun, from dynamo-generated magnetic fields in a differentially rotating late-type star's 
convective zone. For such stars, X-ray luminosity linearly increases with increasing rotation 
until saturation is reached at a relative X-ray luminosity log Lx/L-fo \ ~ —3. The origin of 
this plateau in Lx / -^bcl * s unknown, but possible explanations include saturation of the dy- 
nam o or limits on the num b er and size of ac t ive are as imposed by the surface area of the star 
(e.g. iGiidel k Nazi 120091 ) . iDempsey et al.l ( 1l993aU bl) studied the X-ray emission properties 



of a large sample of rapidly rotating companions in RS CVn systems, where tidal interactions 
amongst the companions helps maintain a more rapid rotation than that of single stars at 
a similar age. The X-ray spectra of the coronae of such stars are best described by 2-T 
thermal plama models comprised of a cool component of a few MK and a hot component 
of tens of MK (IDempsey et al.lll993bl ). EUV spectra of the coronal emission from nearby 
active binary syste ms reveal a continuous em ission measure distributio n (EMD) that varies 
with t emperature ( jSanz-Forcada et al.ll2003l ). The EMD is explained by lCargill fc Klimchuk 
(120061 ) as indicative of a number of unresolved, stable, coronal loops, varying from solar-like 
loops with modest temperatures (~ 2 MK) and densities (~ 10 9 cm -3 ) to hotter loops (> 8 
MK) with higher densities (> 10 9 cm -3 ). 



If we attribute the detected X-ray emission to the main sequence companion in HFG 
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1, then the activity ratio, log Lx/L^^ ~ —3.0 ± 0.2, is at the saturation level. Assuming 
fully ionized coronal plasma, = 0.85n e , the emission measure, EM, determined from 
the X-ray spectral fit is related to the density of coronal plasma by n e = (EM /0.85V) 1 / 2 , 
where V is the emitting volume. The emitting volume can be approximated as a shell on 
the stellar surface by V = 4^i? 3 [(L + l) 3 — 1], where L is the characteristic coronal loop 
length in stellar radii, i?*. For inactive stars, densities up to 10 9 cm -3 and loop lengths 
~ 0.1-R* are typical, while active stars have higher densities, from 10 10 to 10 11 cm" 3 and u p 



to 10 cm in extreme cases, with generally smaller loop lengths, L < 0.1R* (iGiidell 12004 ). 



If we take L = O.li?*, then the lower bound on the density of the hot X-ray gas in HFG 1 is 



~ 10 10 cm 3 , consistent with a corona around the companion. 

Applying similar considerations to the X-ray emission detected from DS 1, if we attribute 
the X-ray emission to the companion, the resulting value of log Lx/L-^^ ~ —2.0 ± 0.3 is 
significantly above saturation. However, the uncertainty in Lx/L^ ^ does not include the 
poorly constrained uncertainties in the temperature and radius of the companion, which can 
raise L-^ \ enough to bring Lx/L-^^ down to the saturation level. Nevertheless, a value of 
Lx/£ DO i as large as ~ —2.0 would not be without precedent. Indeed, many studies find 
a large scatter in t he Ly/L^^ ratio in the saturated region (e.g., the pre-main sequence 
sample studied by IStassun et al.l 12004 ). with some stars reaching log L x /L^ ^ ~ —2.0. 
The possibility of X-ray flaring activity from DS 1 may also account for the high value of 
Lx/L-fo Q \. Applying similar coronal assumptions to the X-ray emission detected from DS 1 
results in a lower bound on the coronal density of ~ 10 10 cm -3 . 

Although the variability of LoTr 5 (in light of the uncertainty in its orbital period) 
suggests that the X-ray emission detected from LoTr 5 might arise from accretion processes, 
there is corroborating evidence that supports the notion that the X-rays arise as a con- 
sequence of coronal emission from the giant companion. Good indicators of such coronal 
activity are the strong Ca II H & K, Mg II, and Ha emission lines, which are likely due to 
chromospheric activity associa ted with the giant com panion than any accretion-related asso- 
ciated with the sdO primary ( jjasniewicz et al.lll996l ). Indeed, the ratio log Lx /L-^^ ~ — 5 
and projected rotational velocity v sin i ~ 67 km s -1 ( jStrassmeier et al.lll997l ) of the giant 
companion in LoTr 5 are consistent with those of rapidly rotat ing, intermedia te mass, G- 
and K-ty pe giants with coronal activity, both in single systems (jGondoinl 120051 ) and binary 
systems (jGondoinl 120071 ) . Suc h giants do not follow t he empirical activity-rotation relation, 
Lx ~ 10 27 (v sin i) 2 , found by Pallavicini et aL dl98ll) for m ai n sequenc e late- type stars and 
do not appear to reach a saturation level ( jGondoinl 120071 ). iGondoinl ( 120071 ) interpret the 
linear relation between the coronal radiative flux density and the average surface magnetic 
flux density as evidence for increasing magnetic surface flux in response to an increasing 
angular rotation velocity. The X-ray properties of LoTr 5 also agree with those by found 
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by iDempsey et al.l ( I1993al lbl) for RS CVn systems. The available evidence indicates that the 



binary in LoTr 5 may be similar to such systems, as also suggested by IStrassmeier et al. 
f )l997h . 



5. Summary 

The X-ray emission sources detected at the binary central stars of three PNe — two 
suspected post-common envelope binaries, HFG 1 and DS 1, and the binary LoTr 5, whose 
period is unknown — reveals properties consistent with coronae in late-type, spun up compan- 
ions. For HFG 1 and DS 1, the measured ratio of X-ray to bolometric luminosity, Lx/L-^ Q ^, 
indicates that the main sequence companions are at or beyond the saturation level deter- 
mined for rapidly rotating late-type stars. The X-ray light curves suggest a steady source in 
HFG 1 and a possibly flaring source in DS 1, which may account for the large relative X-ray 
luminosity measured for the latter central star, if the companion is the source. The value 
of Lx/L-fo Y inferred for the companion star in LoTr 5 is similar to that of rapidly-rotating 
giants and RS CVn systems. The X-ray temperatures and emission measures determined 
from the X-ray spectral fits are furthermore consistent with those found for active, late-type 
stars. Although, the X-ray temperatures would also be consistent with those expected from 
accretion onto the pre- WD, sdO, primary stars, the lack of evidence for Roche lobe overflow 
from the late-type stars in HFG 1 and DS 1 casts doubt on such an origin for the X-ray 
emission. Since the period of the binary in LoTr 5 is unknown, we cannot determine if the 
late-type giant secondary in LoTr 5 is capable of filling its Roche lobe, and Roche lobe over- 
flow remains a possible source of accreting material in this case. However, there is compelling 
and corroborating evidence that the X-ray emission is due to coronal activity associated with 
the rapidly rotating giant. We conclude that the observed X-ray emission in each of these 
three binary CSPN systems is most likely due to coronal activity associated with the late- 
type companions. Although it is clear that supporting observations are required to detect 
and characterize any putative companions more directly, these X-ray observations thereby 
demonstrate the potential utility of X-ray searches for faint companions in the central stars 
of PNe. 
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Facilities: XMM (EPIC), CXO (ACIS). 
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Table 1: Summary of Binary Properties' 



PN 


Mi,M 2 


-Ri, -R2 


T eff,i' T eff,2 


L bol,i' L bol,2 


References 




(M ) 


(^0) 


(kK) 


(A=>) 




HFG 1 


0.57,1.09 


0.19,1.30 


83,5.4 


1500,1.3 


1 


HFG 1 


0.63,0.41 


...,1.15 


...,5.3 


...,0.9 


2 


DS 1 


0.63,0.23 


0.157,0.402 


77,3.4 


776,0.02 


3 


LoTr 5 


0.6,1.1 


0.05,8-12 


185,5.3 


2600,100 


4,5,6 



a The stellar properties of the primary and secondary components are labeled by the 1 and 2 subscripts, 
respectively. 

References: 1-Shimanski et al. 2004, 2-Exter et al. (2003) solution for K p = 49 km s" 1 and i = 29°, 
3-Hilditch et al. (1996), 4-Feibelman & Kaler (1983), 5-Jasniewicz et al. (1996), 6-Graham et al. (2004) 



Table 2: Summary of X-ray Observations 



Object 


Date 


ObsID 


Chip 


^exp 


Net CR 










(ks) 


(cnt ks -1 ) 


HFG 1 


2008 Dec 11 


9954 


ACIS-S3 


11.3 


12.4 


DS 1 


2009 Jul 19 


9953 


ACIS-S3 


23.8 


2.30 


LoTr 5 


2002 Jun 06 


12850201 


EPIC pn 


20.9 


34.7 








EPIC MOS1 


25.7 


11.3 








EPIC MOS2 


25.7 


8.78 


LoTr 5 


2002 Dec 04 


3212 


ACIS-S4 


27.7 


26.1 



Table 3: Spectral Fitting of the X-ray Emission from Binary CSPNe 



PNe x 2 (d.o.f. a ) log N H kT T x log norm. 6 log EM b 


lo § F x,ohs 


log -P^.unabs 


log L x 






(cm 2 ) 


(kcV) 


(MK) 


(cm 5 ) 


(cm 3 ) 


(erg cm -2 s _1 ) 


(erg cm~ 2 s _1 ) 


(erg s _1 ) 


HFG1 


0.79 (510) 


21.4 


5 +,u 
u -°-o.i 


5.8 


4 vH>.l 


53.10 


-13.57 


-13.13 


30.51 










21.3 


4 ^+°- 2 

4 -°-0.2 


53.11 


-13.56 


-13.39 


30.25 


DS 1 


0.49 (510) 


20.9 


o.3l8:? 


3.1 


r 7+0.2 


52.07 


-14.65 


-14.45 


29.32 








1 O+0.4 
x - z -0.2 


14.5 


r o+0.2 
-O.O_ 2 


52.47 


-14.14 


-14.06 


29.71 


LoTr 5 


1.34 (34) 


0.00 


U - DO -0.05 


7.6 


_4 OQ+0.05 


52.54 


-13.54 


-13.54 


29.93 


(XMM) 






9 97+0.51 
z - z ' -0.35 


26.3 


07+O.O6 
"*'°'-0.05 


53.11 


-13.22 


-13.22 


30.25 


LoTr 5 
(CXO) 


0.50 (23) 


0.00 


64+ 07 

u - u ^-0.07 
o 4 q+0.83 
°-^ y -0.64 


7.5 
40.5 


_4 71+0.10 

1 -0.09 
o qo+0.04 
°- yo -0.05 


52.77 
53.55 


-13.30 
-12.74 


-13.30 
-12.74 


30.18 
30.73 



"When fitting the unbinncd spectra of HFG 1 and DS 1, the empty bins are used as constraints, hence the 
inflated number of degrees of freedom (d.o.f.). 

fc Emission Measure (EM) is determined from the normalization value of the mekal model according to: 



EM 



/ 



n e n H dV = 1.2 x 10 58 cm 



norm, x 



D 



1 kpc 



(4) 



-23 - 



0.020 



0.015 

> 

CD 

1/3 
1/3 

e 0.010 h 

o 
o 



0.005 



0.000 
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Fig. 1. — X-ray spectrum of HFG 1 observed by CXO overlaid with the best-fit model 
(top) and residuals of the spectral modeling (bottom). The best fit model (dashed line) 
requires two thermal plasma components at 6 MK and 21 MK, with intervening absorption 
of N H ~ 2.4 x 10 21 cm" 2 . 
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Fig. 2. — X-ray spectrum of DS 1 observed by CXO overlaid with the best-fit model (top) 
and residuals of the spectral modeling (bottom). The best fit model (dashed line) requires 
two thermal plasma components at 3 MK and 15 MK, with intervening absorption of Nh ~ 
8 x 10 20 cm" 2 . 
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Fig. 3. — X-ray spectra of LoTr 5 observed by the three EPIC detectors (pn - circles, MOS1 
- triangles, MOS2 - squares) overlaid with the best-fit model (top) and residuals of the 
simultaneous spectral modeling (bottom). The best fit model (broken lines) requires two 
thermal plasma components at 7.6 MK and 26 MK. No intervening absorption component 
is required, reflecting the small distance and high galactic latitude to LoTr 5. 
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Fig. 4. — X-ray spectrum of LoTr 5 observed by CXO overlaid with the best-fit model (top) 
and residuals of the spectral modeling (bottom). The best fit model (dashed line) requires 
two thermal plasma components at 7.5 MK and 40 MK, with no intervening absorption. 
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Fig. 5. — Unfolded X-ray spectra of LoTr 5 observed by XMM EPIC pn (triangles) and CXO 
ACIS-S (circles) overlaid with the best-fit, two temperature plasma models (CXO - dashed 
line, XMM - dotted line). An increase in flux during the six month interval between the 
XMM and CXO observations is apparent across the energy range. 
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Fig. 6. — Background-subtracted X-ray light curves of HFG 1 (upper panel) and DS 1 (lower 
panel). The light curve of HFG 1 is consistent with a steady count rate, while DS 1 exhibits 
a marginally significant increase in the count rate during the second half of the observation. 
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Fig. 7. — X-ray light curves of LoTr 5 from XMM EPIC pn (upper panel) and the CXO 
ACIS-S (lower panel) observations. 



